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MUTUA'’S METHOD: PARTICULAR POLYNOMIAL AND COEFFICIENTS OF 
PARTICULAR INTEGRALS OF THE FORM ax"e* 


SAMUEL FELIX MUTUA 


Abstract. In many problems involving solutions to ordinary differential equations, students and 
researchers in the field of applied mathematics are faced with challenges of obtaining the particular 
integrals. In particular, when the defined problem is formulated with non-homogeneous Ordinary 
Differential Equations with constant coefficients, Methods of Undetermined coefficients prove to be 
lengthy. Solutions to this will be reached by use of a polynomial, which we shall call particular polynomial 
whose values at given instants will be used to determine the coefficients of these particular integrals. We 


will define the polynomial p,(t) of degree q for q <n ,n been the order of the given ODE, by 


p,(@)=(@-—m,\t-m,)..C t-m,,t-m,).m, being the root of the characteristic function all not 


q-l 
equal to c (in the exponential). If k are the number of roots all equal to c , then the particular integral is, 


cx 


SE k=01,2,a n 


1.0 INTRODUCTION 


Previously the methods of undetermined coefficients and variation of parameters have 
been used as techniques for solving these particular integrals. The former involves 
differentiating the assumed particular integral, substituting to the given ODE and then 
equating the LHS and RHS coefficients to solve for values of the coefficients. In the latter 
method, the constants are replaced with parameters (functions) which are varied using 
the given ODE to solve them. Use of the particular polynomial will simplify this process 
of obtaining such particular integrals. This will prove helpful as it will give results 
without having to go through the entire complex algorithms encountered in the past. 


2.0 MUTUA’S METHOD 
2.1 PRELIMINARIES 


Consider the second order ODE 


2 
dy d 
AE EE ae 
dx ` dx 





Its characteristic function is 
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m°’ +3m+2=0 


Which has roots m, = -1 and m, = -2 and the complementary function is 


To solve for the particular integral using method of undetermined coefficients, since 2 is 


nota root we let 
Vane ,y pee ae yy So Ae” 


Substituting back to the given ODE yields 


(2°) Ae ** + 3(2') Ae “* + 2Ae ** = 12 Ae ** = e° (RHS) 


Now if we critically look at the coefficient — , its denominator has a very close 
12 


correlation with the value of characteristic function f(m) = m* +3m +2 atm =2 


f(2)=27 +3(2)+2=12 


2.2 THE PARTICULAR POLYNOMIAL 


Now, we consider a general linear ODE of order n with constant coefficients, 





d"y d” y dy z 
a, +a, , tg Gt OY ae (1) 
dx” dx" dx 


Where a, « 0 and cis a constant, 
We have the characteristic function 


n-l 


n 
f(m)=a,m +a, m +.. +a,m+a, (2) 


We define a polynomial p, (t) of degree q for q < n by 


p,(t)=(t-m,\t-m,).( t-m,,Xt-—m,) (3) 


MUTUA’S METHOD 
Where m,,m,,..., m, is a subset of the roots of f (m) alls c. 


We shall call this particular polynomial as it will be used to determine the particular 
integral. 


CASE1: If all the roots of the characteristic function are not equal toc, we find that 
f(c) +0 


And therefore we take our particular polynomial to be 














p(t)=a,m"+a_.m" +..+am+a, (4) 
Now, the coefficient 4 of the particular integral, y, = Ae“ becomes A = i.e. to say 
p(c) 
. : 1 
that Particular integral y, = e“ CASE2: 
p.c) (5) 
If m, m, a m, c and M, ọso M, #C 
Then we realise that; 
f(m) =(m=m, Ym —m,).< m= mom =m, m—=m,) (6) 


Yields f(c) = 0 „Since the first k factors become zero. 


1 
f(c) 


In this case the coefficient of particular integral cannot be 





. Clearly it can be seen 


that e“ will be absorbed in the complementary function. 


eX 2- vex k-1 cx 


Also, xe“, x°e",.., x° e° are contained in the complementary and they cannot form the 


particular integral. 


cx 


In this case we assume a particular integral of the form y, = Ax “e 


(7) 
Remark 


k cx 


From the differentiation of x*e™ (by product rule) where k e z` we realise that as the 
order of differentiation increases, in some of the terms the value of finite k diminishes 


and ultimately tends to zero. The first term with k =0 is in the k” ordered derivative 


x 


and will be of the form Be “ where 


B=k! 


J 


It must be noted that if the method of undetermined coefficient was to be used the 


th 


resulting k” ordered derivative needs to be differentiated (n - k) more times so that 
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upon substitution to the given ODE and equating the LHS and RHS coefficients 
particular integral is obtained. 


Now; 

When k =1, e“ is absorbed in the complementary function. 

Wetake y, = Axe “ and the particular polynomial as 

p(m) = (m—m,\m—m,)..( m-m,) (8) 
Where p(c) # 0 


Using the knowledge of integration by parts (its reverse), we arrive at the following. 





The coefficient A as 
1.p(c) 


When k = 2, e“, xe“ are absorbed in the complementary function. 


cx 


We take y,= Axe 


P(m)=(m—-—m,)\m—m,)..0 m-—m,) (9) 


Once again p(c) # 0 


1 1 


The coefficient 4 becomes A = ———— = 
2.1.p(c)  2!p(c) 





For k = 3,4, q <n 


1 1 





PEREN 


The corresponding coefficients of A are; : 
BIP te) AP, key. aP, 


x 4 cx pm: 


And the particular integrals are Ax *e”, Ax *e™.... Ax “e”, with A, as defined above. In 


general, the particular integral is given by; 


1 l 
i a a E EAN, (10) 
k! P, (c) 


Equation (10) illustrates MUTUAS METHOD of finding the particular integral. 


Note: If the RHS of equation (1) is multiplied by a constant, then equation (10) is also 
multiplied by the same constant when evaluating the particular integral. 


MUTUA’S METHOD 
3.0 APPLICATION 


The above method can be used in modelling phenomenon in which the governing 
mathematical equations are linear non-homogeneous Ordinary Differential Equation 
with constant coefficients, the RHS being an exponential function. 


3.1 NUMERICAL TEST CASES 


EXAMPLE 1 





d 
Solve 2——+5 3y =e. 
dx dx 


Solution 

The characteristic function is 2m° + 5m -3 = 2(m +3Xm -1/2)=0 
m, = -3,m, =1/2 

.. Complementary function y, = Ae ** + Be 
Here all roots are not equal to 3 


.. Particular polynomial P (m) = 2m* + 5m —3, P(3) = 2(3)'5(3) -3 = 30 


3y 


1 
Thus E 
30 


And the complete general solution of the given ODE is 


EXAMPLE 2 


d° d’°y dy 
Solve 7 Lag aia eee 
dx dx dx 


x 





The characteristic function is m° - m° -4m +4=(m-1Xm-2Xm+2)=0 








Here m, =2=c¢ 


.. Particular polynomial P (m) = (m -1)X(m + 2), P(2) = (2 -1X2 +2) =4 
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EXAMPLE 3 


d° d’ dy s 
Solve aag ee y = 4e” 
dx” dx dx 





The characteristic function is m° —3m* +3m-—1=(m-1)* =0 











.. Complementary function y, = Ae * + Bre * + Cx *e* 





Here m, m, m, l=c 





.. Particular polynomial p(m) =1, P(1) =1 


E E 2 


4 2 
Thus y= — xe =—x’e 
3 


And the complete general solution of the given ODE is 


2 3 





y =y, +y, = Ae" + Bxe*+Cxr’e*"+5x e` 
3 
EXAMPLE 4 
d’ y dy 
Solve —+2—+5y=se" 
dx dx 


Solution 


The characteristic function is m? +2m +5 =(m+1+2im+1-2i) =0 
m = -1-2i,m, =-14+2i 
.. Complementary function y, = Acos x + B sin x 


Here all roots are not equal to -7 


<. Particular polynomial P(m) = m° +2m +5, P(-7) = 2(-7)° + 2(-7)+5 = 40 


MUTUA’S METHOD 


-7x 


1 
Thus WS w 


And the complete general solution of the given ODE is 


1 
y= y., +y, =Acos x+ Bsin age 


Tx 


Note: The above procedure can also be used when c ¢ Z 


EXAMPLE 5 





d’ 
Solve Z554 6y= 2e" 


Solution 


The characteristic function is m* -5m +6 =(m-3Xm-2)=0 





m =2,m, =3 
.. Complementary function y, = Ae °" + Be ` 
Here all roots are not equal to -1/3 


.. Particular polynomial P(m) = m° -5m +6, P(-1/3) = (-1/3)° —5(-1/3) +6 = 70/9 





-1/3x 9 -1/3x 
=e 


9 
Thus y, 2 Je = 
70 


COMPARISON OF MUTUA’S METHOD OF PARTICULAR POLYNOMIAL WITH THE 
METHOD OF UNDETERMINED COEFFICIENTS 


Suppose we redo example 3 with the METHOD OF UNDETERMINED COEFFICIENTS 


a” d’ dy , 
Solve = ie ee y = 4e" 
dx” dx dx 





The characteristic function is m° -3m° +3m-1=(m-1)* =0 
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.. Complementary function y, = Ae * + Bre * + Cx *e* 





To solve for the particular integral, we assume a particular integral of the form. 


Where 4 is a constant to be determined. Differentiating the above equation three times 
we get, 


y's =3Ax e” + Ax e” 
y, = 6 Axe ` +3Ax e” +3Ax e” + Ax e” =6Axe” +6Ax e` + Ax`e” 


y''", = 6Ae* +6Axe ` +12 Axe * +6Ax e" +3Ax°e* + Axe" =6Ae* +18 Axe * +9Axe* + Axe” 
Substituting into the given ODE with have 


(6 Ae * +18 Axe * +9 Axe" + Ax “e*)—3(6Axe * + 6Ax e” + Axe") +3(3 Ax e” + Ax?) Ax e” = 4e” 
Or 


6Ae = 4e 
2 
A=— 
3 
2. ge 
>y,= xe 
3 


And the complete general solution of the given ODE is 


3.2 CONCLUSION 


Clearly from the above example we conclude that MUTUA’s method of particular 
polynomial arrives to the required particular integral of the form Ax"e* faster than the 
already existing alternatives. In particular when the ODE in question is of high order a 
lot of time is spent on differentiation. Thus MUTUA’S method is of great significance 
since it simplifies the computation process involved in obtaining particular integrals. 
This is an important aspect to any mathematical formula. 
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